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ABSTRACT 

Recently Gizon, Duvall, and Schou (2002) suggested that supergranulation 
has a wave-like component. In this paper I show that the same phenomenon can 
be observed using surface Doppler shift data, thereby confirming their observa- 
tions. I am also able to measure the dispersion relation to lower wavenumbers 
and to extend the results for rotation and meridional flows beyond ±70° latitude. 



1. Introduction 

Solar supergranulation, which was first described more than 40 years ago by Leighton, 
Noyes, and Simon (1961), has remained difficult to explain. While it appears to be similar to 
the granulation, several observations, such as the anomalous flows obtained using correlation 
tracking or similar methods (e.g. Beck and Schou (2000), in the following BS), challenge the 
interpretation as a simple convective phenomenon. An alternative explanation is that the 
supergranulation has a wave-like component and observational evidence for this was recently 
reported by Gizon, Duvall, and Schou (2002) (GDS in the following). 

It has been the cause of some concern that the wave like character has not been observed 
previously, given the many measurements of the rotation rate of the supergranulation and 
attempts to determine why it appears to be anomalous. A likely explanation is that it was 
obscured by the limitations of the surface Doppler shift data. 

Supergranular motions are predominantly horizontal near the photosphere. Let be 
the velocity in the longitudinal (0) direction and Vq the velocity in the latitude (6*) direction. 
The observed Doppler velocity V is then given by 

V = Vtf) sin (f) + Ve cos (f) sin 9 (1) 

under the assumption that the observations are made in the equatorial plane (i.e. that BO 
= 0) at a large distance from the Sun. For the component, in particular, the projection 
leads to a large modulation of the signal as the Sun rotates and to a significant spreading of 
the power in the frequency domain. 
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2. Data Analysis and Results 

The data used here are the same as those used in BS, which provides further details. 
They consist of a 60 day sequence of 1024^ full disk Doppler velocity images obtained by the 
MDI instrument (Scherrer, et al. 1995) on the SOHO spacecraft, covering the time period 
1996 May 24 to 1996 July 22. These images were derotated and averaged to a 60 minute 
cadence, remapped to a uniform grid in longitude and latitude and binned down to have a 
resolution of 1200 points covering 360° of longitude and 600 points in latitude between ±90°. 
A temporal average was removed as part of this process. 

Strips in longitude were extracted, multiplied by a weighting function, tracked at a 
latitude dependent rate, zero padded to 360° by 12° and passed through a 3 dimensional 
Fourier transform. The weighting function is a product of a function of 6 and a function 

of (f). The function of ^ is 1 for \59\ < 4.5°, for \59\ > 5.5° and falls off as a cosine bell 
for 4.5° < \d6\ < 5.5°, where 66 is the distance from the target latitude. Two weighting 
functions in ^ were applied. One is optimized for recovering V^, the other for recovering 
Vq. The first weighting function is given by W^p = A^sigrKp/ a/ sin^ + 0.01, where is an 
apodization function which is for |0| = 0° and for |0| > 80°, 1 for 10° < |0| < 70° and 
cosine bell apodized for 0° < |0| < 10° and 70° < |0| < 80°. The second weighting function 
is given by We — Aq/ sjcos^ + 0.01, where Aq is only apodized between 70° and 80°. The 
small bias of 0.01 was added to avoid singularities. 

Figure 1 shows the power as a function of direction of propagation. Note that there is 
significantly more power in the prograde than in the retrograde direction, consistent with 
the findings of GDS. Also there is fittle power in the N-S direction where GDS did observe 
significant power. Since the method used by GDS is sensitive to both components of the 
velocity this indicates that the displacement, at least at the surface, is almost entirely in the 
direction of propagation. 

Given the above waves in the direction should show modulation given by A^ when 
using the weighting and waves in the 6 direction should be modulated by Aq when using 
Wq. For other directions the modulation is more complex and only waves propagating near 
the longitude and latitude directions are studied here. The velocity projection also means 
that waves traveling in latitude are essentially unobservable at the equator. At high latitudes 
the foreshortening becomes a problem. 

Figure 2 shows examples of the resulting power spectra. Note the different power level 
in prograde waves (lower right) and retrograde waves (upper right) in the left panel. In 
the right panel there is excess power in the southbound waves (upper right) relative to the 
northbound waves (lower right). 



-3- 



The advection of a wave pattern will, at a given azimuth, lead to two peaks in the 
power spectra at frequencies u}± = ±uJo{k) + kv, where uJo{k) is the mode frequency, k is the 
wavenumber (given here in terms of the equivalent degree / = kRQ, where Rq is the solar 
radius) and v is the component of the velocity in the direction of propagation. 

To quantify the properties of the waves, power spectra at selected k were extracted and 
a model consisting of two Lorentzians and a background 



was fitted to them using a maximum likelihood algorithm. Here A± are the amplitudes, w 
is the the (common) width of the modes and B is a constant background, The dispersion 
relation may be derived as ujq = (a;+ — ci;_)/2 and the flow velocity from v = (a;+ + u}_)/2k. 
To account for the projection factor and the finite observation time, the model is convolved 
with an estimate of the appropriate window function. 

The power spectra are binned by a factor of 5 in kx, rebinned to a uniform grid in k 
and azimuth and directions within ±15° of the desired direction are averaged before being 
fitted. This averaging ensures that the log of the resulting power spectra is close to normally 
distributed with a uniform variance, allowing for a simple maximum likelihood algorithm 
and stable fits. The inaccuracy in the parameters caused by the averaging over azimuth is 
neghgible. 

Figure 3 shows the dispersion relation for the waves traveling in longitude at the equator. 
The solid hue shows the fit of a power law over the interval 50 < Z < 200, which gives an 
exponent of 0.56. Points below I ^ 30 have proven difficult to fit and may be unrehable. 

Figure 3 also shows the rotation rate at the equator, which is almost independent of 
I, allowing the results to be averaged over I, as done in Figure 4, which also shows the 
meridional flows. 

Finally Figures 5 and 6 shows the amplitude, line width and anisotropy. 



Possibly the main result is to confirm the results of GDS. Given the differing analysis 
techniques this represents a significant confirmation and lays to rest the concerns based on 
the fact that the phenomenon had not been observed previously in surface Doppler shift 
data. 





3. 



Discussion 
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While not physically motivated it is interesting that the frequencies appear to follow a 
square root quite well. As can be seen from Figure 3 the dispersion relation at 40° is quite 
similar to that at the equator in both directions. 

The small variation of the rotation rate with I is consistent with a simple depth inde- 
pendent flow. This is in stark contrast to the numbers derived from the same data in BS, 
which varied by lOnHz between /=50 and /=200, and which were difficult to explain. The 
rotation rate derived here is remarkably close to the magnetic tracer rate of Komm, Howard, 
and Harvey (1993). The similarity may be related to the fact that the advection of the small 
magnetic elements appear to be influenced by the supcrgranulation. This similarity is quite 
unlike the results in BS, in which the rotation rate, at least at low /, exceeded that of the 
plasma at any depth. This lends further credence to the idea that the flows inferred here 
represent physical advection rates. 

The meridional flow results appear to be reliable to at least ±70° latitude with a clear 
turnover between 20° and 30° and no sign of a second cell. Note that the meridional flow 
is not perfectly antisymmetric across the equator. As shown this is largely explained by a 
P angle error caused by the misalignment of the MDI instrument on the SOHO spacecraft 
(Toner 2000) and by the difference between the Carrington elements and the true rotation 
axis of the Sun (Giles 1999). An error in the rotation axis causes part of the rotation velocity 
to be misidentified as meridional flow. 

Figure 5 shows that the power peaks around I — 100, as previously seen. The hnewidth 
appears to be an increasing function of I, possibly indicating that smaller scale features damp 
more quickly. The increase for Z < 50 is likely an artifact of the difficulty in fitting low degree 
features. 

As in GDS Figure 6 shows that the power is highly anisotropic at all latitudes and that 
the anisotropy is in the prograde direction and increasingly towards the equator at higher 
latitudes. 

The physical nature of these waves, if that is indeed what they are, is still unknown. 
However, with the results presented here, it may be possible to further constrain the models. 

Further progress may be possible by extending these results over time or by probing the 
depth structure using time-distance helioseismology. 
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Fig. 1. — Power (averaged over 30 < / < 200 and < v < 10/i.Hz) at the equator as a 
function of azimuth using the Wo weighting. Prograde is at 180° while northward is at 270°. 
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Fig. 2. — Left: Cut in a 3 dimensional power spectrum using PF^ at the equator around 
ky — 0. The grayscale is logarithmic with black 300 times the power of white. Right: Cut 
at 40° latitude around k^; — using Wo- The poor resolution is due to the narrow strips in 
latitude. The power is symmetric with respect to a simultaneous change of sign of I and v. 
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Fig. 3. — Results for waves traveling in longitude at the equator. Left: The dispersion 
relation (plusses), 1.18(V50)° VHz (dotted) and 1.11(//50)°-^VHz (sohd). Also shown are 
the results at 40° for waves traveling in longitude (diamonds) and in latitude (squares). The 
results at high / are less reliable due to the foreshortening. Given the poor resolution in kx 
(see Fig. 2) only the results at every 20 in / arc shown for the waves travelign in latitude. 
Right: The sidereal rotation rate (solid) and the equatorial rate from Komm, Howard, and 
Harvey (1993) (dashed). 
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Fig. 4. — Derived flows averaged over 100 < I < 200. Left: The sidereal rotation rate 
derived using (plusses) and the rotation rate from Komm, Howard, and Harvey (1993) 
(solid). Right: The meridional flow (plusses) (using We) and an estimate of the effect of the 
incorrect P-angle (dotted). 
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Fig. 5. — Power 'w{A+ + A-)/2 and linewidth ty as a function of I for longitudinally travehng 
waves at the equator. 
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Fig. 6. — Anisotropy as a function of latitude averaged over 100 < I < 200. Plusses show 
the ratio of prograde to retrograde power. Diamonds show the ratio of Northbound power 
to Southbound power. The anisotropy is essentially independent of I. 



